Corrugated waveguide under scaling investigation 
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Some scaling properties for classical light ray dynamics inside a periodically corrugated waveguide 
are studied by use of a simplified two-dimensional nonlinear area-preserving map. It is shown that 
the phase space is mixed. The chaotic sea is characterized using scaling arguments revealing critical 
exponents connected by an analytic relationship. The formalism is widely applicable to systems with 
mixed phase space, and especially in studies of the transition from integrability to non-integrability, 
including that in classical billiard problems. 
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Interest in the problem of guiding a light ray inside 
a periodically corrugated boundary has increased in re- 
cent years, partly because the topic is applicable in so 
many different fields of science. Applications involve e.g. 
ray chaos in underwater acoustics [E 0, S], scattering 
of a quantum particle in a rippled waveguide [1], quan- 
tized ballistic conductance in a periodically modulated 
quantum channel [HI] , quantum transport in ballistic cav- 
ities [6], transport through a finite GaAs/AlxGai_xAs 
heterostructure |21], comparisons of classical v. quantum 
behavior in periodic mesoscopic systems jSj, and anoma- 
lous wave transmittance in the stop band of a corrugated 
parallel-plane waveguide [§]. 

There are many different ways of describing problems 
involving waveguides. One of them is the well known 
billiard description, and such a procedure is used in the 
present paper. A billiard problem consists of a system 
in which a point particle moves freely inside a bounded 
region and suffers specular reflections with the bound- 
aries. Generally, the dynamics is described using the 
formalism of discrete maps. Depending on the combi- 
nations of control parameter values, as well as on the 
initial conditions, the phase spaces for such mappings 
fall into three distinct classes, namely: (i) regular; (ii) 
ergodic; or (iii) mixed. Generally, the integrability of the 
regular cases is related to angular momentum conserva- 
tion, the static circular billiard being a typical example. 
On the other hand, for completely ergodic billiards, only 
chaotic and unstable periodic orbits are present in the 
dynamics. Two examples of case (ii) are the Bunimovich 
stadium [l3| and the Sinai billiard [Tlj . For these two sys- 
tems, the time evolution of a single initial condition, for 
the appropriate combinations of control parameters, is 
enough to fill the whole phase space, ergodically. Finally, 
there are many billiards of mixed phase space structure 
[H, d, [3, [H, m, , having a range of control param- 
eters whose physical significances differ. Depending on 
the combination of both initial conditions and control pa- 
rameters, the phase space presents a very rich structure 
containing invariant spanning curves (sometimes known 
as invariant tori), Kolmogorov-Arnold-Moser (KAM) is- 
lands, and chaotic seas. 



In this Letter, I use scaling arguments to describe the 
behavior of the variance of the average reflection angle, 
within the chaotic sea, for a classical light ray undergo- 
ing specular reflections inside a periodically corrugated 
waveguide. The model consists of a classical light ray 
that is specularly reflected between a corrugated surface 
given by y = 2/0 + dcos{kx) and a flat plane surface at 
y = 0. The term yo denotes the average distance between 
the corrugated and flat surfaces, d is the amplitude of the 
corrugation and k is the wave number. The dynamical 
variables used in the description of the problem are the 
angle 9 of the ray's trajectory measured from the posi- 
tive horizontal axis, and the corresponding value of the 
X coordinate at the instant of reflection. Moreover, the 
mapping is iterated when the light hits the surface y = 0; 
thus multiple reflections with the corrugated surface can 
be neglected. 

The problem lies in obtaining the map r(a;„,0„) = 
(xn+i, 0„4_i), given the initial conditions (a;„,0„) shown 
in Fig. [TJa) . From the geometrical considerations illus- 
trated in the figure, it is easy to see that a;* — x„ — 
{yo + dcos{kx'^))/9n, and similarly that 



Xn+i - a;* = (t/o + dcos{kxl))/9. 



n+l 



(1) 



The term gives the location of the collision on the 
corrugated surface. The angle 0„ is written as 



(2) 



where tan('0„(a;)) — dy(x)/dx — ~dksui{kx^) gives the 
slope of the surface at a: = x*. Equations (P) and ^ 
correspond to the exact mapping. However in this paper 
we make the following approximations: (a) we assume 
that d/yo <C 1, so that yo + c?cos(fca;* ) = yo; (b) in the 
same limit, we also assume that tan('0„) = ipn- The 
condition d/yo 1 holds in practice for studies of wave 
propagation in a corrugated waveguide [3] . It also applies 
to the investigation of transport in mesoscopic channels 
8] and can even be used in the characterization of wave 
propagation in deep ocean P, 0] ■ In the description of the 
mapping, it is convenient to use dimensionless variables, 
6 = d/yo, "fn = On/k and Xn = kxn/yo. Thus, the 
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FIG. 1: (a) Reflection from the corrugated surface of a light 
ray coming from the flat surface at i/ = 0. (b) Details of 
the trajectory before and after a collision with the corrugated 
surface. 



simplified two-dimensional nonlinear mapping is given by 



T 



Xn+l — Xn 



1 

7>i + i 



mod(27r) 



7„+i = 7„ + 2S sin X 



(3) 



The determinant of the Jacobian matrix for this map- 
ping is unity, so that it is area-preserving. The mapping 
([3]) is described by a single effective control parameter, 
d. US — 0, the system is integrable while for <5 / 0, the 
phase space is mixed, containing both chaos and regular- 
ity (fixed points and quasiperiodic behavior). The system 
thus experiences an abrupt transition from integrability 
to non-integrability when going from S — to S ^ 0. 
I shall investigate some dynamical properties for very 
small values of S (which meets the initial approximation 
d/yo ^ 1) near the transition, using scaling arguments. 
There is a wide range of systems exhibiting transitions 
from integrability to non-integrability, and the formalism 
used in the present model can be in principle be extended 
to them, in particular to identify classes of universality. 

Fig. [IJa) shows the phase space generated by iterating 
the map ^ for the control parameter 5 = 10~^. It is 
easy to see that the phase space is of mixed form, includ- 
ing invariant spanning curves, KAM islands and a large 
chaotic sea for both positive and negative values of the 
variable 7. The invariant spanning curves limit the size of 
the chaotic sea and correspondingly yield a limited range 
of possible values for the chaotic orbits. Fig. [D(b) shows 
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FIG. 2: (Color online) (a) Phase space for the mapping (|3]) 
with the control parameter 5 — 10"''. (b) Positive Lyapunov 
exponent obtained via the triangularization algorithm. 



positive Lyapunov exponents obtained for the chaotic sea 
using the triangularization algorithm |18l] , for the control 
parameter 6 = 10~^. An ensemble of 10 different initial 
conditions was iterated 5 x 10* times. The initial condi- 
tions considered were Xq ~ S and 10 different values of 
70 uniformly distributed in the interval 7 e [0, 27r). The 
average value obtained was A — 1.624(7), where 0.007 is 
the standard deviation obtained from the 10 samples. 

The existence of invariant spanning curves limiting the 
size of the chaotic sea confers an interesting property on 
the chaotic orbits. From an initial condition in the re- 
gion of the chaotic sea, the system wanders within the 
entire accessible region but it is always confined between 
two invariant spanning curves. The location of the first 
positive and negative invariant spanning curves depends 
on the value of the control parameter. As a consequence, 
the "amplitude" of a chaotic time series is also depen- 
dent on the control parameter, so that the average value 
tends only gradually towards a regime of convergence. 
With this property in mind, I now explore the behavior 
of the deviation of the average value for the angle 7, or 
roughness, an observable defined as 



1 *^ 

uj(n, 5) = — > 



7,2 (n, (5) - 7^2(^,5) 



(4) 
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where 

1 " 

lin,S) = -J2j.. (5) 

i=0 

Equation ([4]) was iterated using an ensemble of M = 
5 X 10^ different initial conditions. The variable 70 was 
kept fixed at 70 = 10~^d while the 5000 values of 
were uniformely distributed along Xq € [0, 27r). Fig.[3I^a) 
shows the behavior of three different values of the rough- 
ness for different control parameters. One can observe 




FIG. 3: (Color online) (a) Roughness a; as a function of the 
iteration number n for three different control parameters, (b) 
After a suitable change of variables, all three curves start by 
growing together, (c) Their collapse onto a single universal 
plot. 

that, as the number of iterations increases, the roughness 
grows at first, but then suddenly bends over towards a 
regime of saturation. The changeover from growth to sat- 
uration is marked by a characteristic crossover iteration 
number denoted as rix- Note that different values of the 
control parameter generate different roughness curves for 
small iteration numbers n. Applying the transformation 
n nS^ yields a coalescence of all the curves at small 
iteration number, as can be seen in Fig. ^h). Based 
on the behavior of Fig. [3]Ja), the following three scal- 
ing hypotheses are proposed: (i) that for small iteration 
numbers, say n ^ rix, the roughness behaves according 
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FIG. 4: (Color online) Plots of (a) i^sat and (b) the crossover 
iteration number 

Tlx SI'S function of the control parameter S. 
to the power law 

uinS^J) (x{nSy , (6) 

where we refer to the critical exponent /? as the growth 
exponent; (ii) for large enough iteration number, Ux, 
the roughness approaches a regime of saturation marked 
by a constant "plateau" given by 

Wsat OC (5" , (7) 

where we call the exponent a the roughening exponent; 
finally (iii) the number of iterations that characterizes 
the crossover, i.e. that marks the change from growth to 
saturation is written as 

Ux (X S"" . (8) 

The exponent 2; is a dynamical exponent. These three 
scaling hypothesis are extensions of the formalism used in 
surface science (see e.g. Ref. 19]). Based on these three 
initial suppositions, the roughness can now be described 
in terms of a scaling function of the type 

Lu{n6'^,S) ^lij{l''72S^,l''S) , (9) 

where / is a scaling factor, and a and b are scaling expo- 
nents. Moreover, the exponents a and b must be related 
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to the critical exponents a, (3 and z. Because I is a scaling 
factor, we can specify that I = (n(5^)'^~^/°^ and rewrite 
Eq. jH as 

uj{n6'') = (n(52)(-i/a)^^([„52j(-b/a)^) ^ (^q^ 

where the function lji ~ a>(l, [n(5^]^''/''(5) is assumed to 
be constant in the limit oi n <^ n,j.. Comparison of Eqs. 

and (O allows one to conclude immediately that 
13 = -l/a. Choosing now that I = 5~^/^, Eq. dH) is 
given by 

oj{n5^,5) = r^'''oj2{5-''''/^^n5^) , (11) 

where the function L02 = Lu{S^^"'''''^nS'^ , 1) is supposed 
constant for n ^ rix- Comparing Eqs. (fTTj) and ([7]), we 
find that a = —1/6. Given the two different expressions 
for the scaling factor one must conclude that the rela- 
tion between the critical exponents takes the form 

. = ^-2. (12) 

Note that the scaling exponents are all determined if 
the critical exponents a and /3 can be obtained numer- 
ically. Fig. [3)Ja) and (b) show the behavior of uJsat and 
rix respectively as functions of S. The saturation val- 
ues were obtained by extrapolation because, even after 
almost lO^rix iterations, the roughness is still not ap- 
proaching its saturation value. Fitting a power law to 
the results of Figs. HUa) and (b), it was found that 
a = 0.509(1), f3 = 0.4997(8) and z = -0.979(4). The dy- 
namical exponent z can also be obtained from Eq. (|12p . 
by evaluation of the numerical values of a and /3, yield- 
ing z = —0.9814(3). This result is in excellent agreement 
with the numerical result obtained in Fig. Hfb). 



With the values of the critical exponents now obtained, 
the scaling hypotheses can be verified. Fig. [3]Jc) illus- 
trates the collapse of the three different roughness curves 
generated from different values of the control parameter 
onto a single universal roughness plot. Such a collapse 
confirms that the scaling suppositions are indeed correct 
despite the three decades used in the control parameter S. 
In the construction of all the roughness curves, the value 
of the dynamical variable 70 was set at 70 = 10^^(5. If 
one increases the initial value of 70, the scaling hypothe- 
ses will remain valid, but an extra scaling time (say n'^) 
appears. Even considering the new extra crossover it- 
eration number, scaling exponents can still be obtained 

In summary, the problem of light ray reflection in a 
periodically corrugated waveguide has been addressed. 
The chaotic sea confined between the first two invari- 
ant spanning curves (the positive and negative ones) was 
characterized using scaling arguments. In particular, the 
deviation of the average value of 7 was obtained as func- 
tions of both the control parameter S and the iteration 
number n. Three scaling hypotheses were proposed, and 
then confirmed by the perfect collapse of all curves onto a 
single universal plot, confirming that the system experi- 
ences a transition from integrability to non-integrability. 
Finally, the obtained critical exponents a, (3 and z allow 
us to conclude that the present corrugated waveguide 
model belongs to the same class of universality of the 
one-dimensional bouncing ball model. 
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